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Abstract 

^ ' We treat 0+ Oi^/3/3 decays taking into account recoil corrections to the nuclear 

■ currents. The decay probability can be written as a quadratic form of the effective 
[ coupling constants of the right-handed leptonic currents and the effective neutrino 

Q\ I mass. We calculate the nuclear matrix elements for the 0"*" — > 2f Qvf3f3 decays of 

■ '''^Ge and ^'^''Mo, and demonstrate that the relative sensitivities of 0"*" 2+ decays 
, to the neutrino mass and the right-handed currents are comparable to those of 

^ \ 0+ — > 0"'" decays. 

ix. 

^ . The neutrinoless double beta (Oz//5/5) decay can take place through an exchange 

of neutrino between two quarks in nulei if the electron neutrino is a Majorana 
\ particle and has a nonvanishing mass and/or right-handed couplings [1-3]. 

^ ' There may be other possible mechanisms such as those involving supersym- 

metric particles which also cause the decay of two neutrons into two protons 
and two electrons [4-6]. In the present work, however, we restrict ourselves 
to the conventional two-nucleon and A mechanisms of Oz//?/3 decay through 
light Majorana neutrino exchange. From the analises of experimental data on 
0"*" — ^ 0^ 0z//9/3 decays, stringent limits on the effective neutrino mass and the 
effective coupling constants of the right-handed leptonic currents have been 
deduced (see e.g. [3,7] and the references quoted therein). On the other hand it 
still seems to be believed widely that 0^ 2+ OujSjS decays are sensitive only 
to the right-handed currents. In view of the theorem that the electron neu- 
trino should have a nonvanishing Majorana mass if Ou/SlS decay occurs anyway 
[8-10], an observation of Oz//3/5 decay due to right-handed interactions would 
certainly mean also a nonvanishing Majorana mass of the electron neutrino. 
The purpose of the present work is, however, not to investigate the role of the 
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Majorana neutrino mass in such a sense, but to demonstrate that it causes 
0+ 2+ Qupp decays directly. 



A direct contribution of the neutrino mass to 0+ — > 2+ 01/(3(3 decays was 
considered in [3] taking into account the nuclear recoil currents, and the inverse 
half-life was given as 

[ri%(0+ - 2+)]-^ = + Fi-(^i-)' + F2+(Z,^r + F24Z2-r , 

(1) 

where Fj± (j = 1, 2) are the phase space integrals and 
Zi± = M^{X) - M^{ri) ± , 

IJle 

Z2± = M:^{rj)±M^^, (2) 
me 

with the electron mass me and 



j 

{v)-vT!UejVej. (3) 
j 

Here rrij is the mass of the eigenstate Majorana neutrino Nj. Uej and Vej are 
the amplitudes of Nj in the left- and right-handed electron neutrinos, A and 
r] the coupling constants of the right-handed leptonic current with the right- 
and left-handed hadronic currents, and the summations should be taken over 
light neutrinos {rrij <S 100 MeV). The nuclear matrix elements {a — A, 77, 
m) are defined by 

= m\Y.^:r:^i^Xmpt) ■ (4) 



The explicit forms of the two body operators M\, and M'^ were given in 
[11] including the contribution of the A mechanism, in which the Ou(3(3 decay 
proceed through an exchange of a Majorana neutrino between two quarks in 
the same baryon in a nucleus. On the other hand the operator iVfm was derived 
in [3] as 



(-^m)nm — ~|wic|[?"nm ® {<^nCm — <T^C„)] 



(2) 
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+i{gw/gA)[rnm ® {Dn X (Tm- DmX <T„)]^^) 

+ {9y/9A?Km ® {D^ - Dm)f^}H{rnm) , (5) 

where = 'f'n — f'm-i H{r) is the neutrino propagation function, and 
gp^ the vector and axial vector couphng constants. C„ and £>„ are the recoil 
correction terms to the axial vector and vector nuclear currents [2,12] given 

by 



Cn = (Pn+Pn)-^^n/2M, 



(6) 



where and p'^ are the initial and final nucleon momenta, M the nucleon 
mass, and ixp — 4.7. The above expression for is, however, not suitable 
for numerical calculations as it stands. Therefore, as was done for Mx, iW^ 
and M'^ in [11], we expand it in terms of the operators Minm with simpler 
spin and orbital structures. 



(7) 



We define the matrix element Mi of the operator Minm analogously to Eq. 
(4). The coefficients Cmi and the two-body operators Minm are listed in Table 
1, where 



h — fnmH(Tnm) i 
Sxnm = [o-n® O-m]^^^ , 

y Knm [^nm ® '''nm] i 

I'+nm — ~l~ I'm : 
Pnm 



Knrn [^nm ® '''+nm] ^ ^ ij'+nm/ ^nm^ i 



Knm, 



iWr. 



]{K) 



a — a/\a\ , 

Pnm — Pn "I" Pm ■ 



(8) 



As was described in detail in [11], Mx and Ai"^ can be expanded in terms of 
Minm with 1 < i < 5, 8 < i < 13, and M'^ in terms of Ai"j„mWith i = 6, 7 
(for the definition of Minm with 6 < i < 13, which do not appear in Table 1, 
see [11]). Of these operators, Minm with 8 < i < 13 are related to the Oz//9/3 
transitions which involve virtual A particles in nuclei, and they are induced 
by the operator M2nm interpreted as acting on two quarks in a nucleon or a 
A particle. 

The new operators Minm with 14 < i < 25 appear only in the expansion of 
Mia- In the derivation of Cmi hsted in Table 1, we have not taken into ac- 
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count the A mechanism yet. Under the same assumption of the non-relativistic 
constituent quark model about the A mechanism as was made in [11], the op- 
erators in Table 1 except Minm with i = 2, 16, 17 do not contribute when 
interpreted as acting on two quarks in a nucleon or a A particle. Since the 
relation fxp — qy — — ^ holds for the quark currents, we see Cmj = for 
i = 2, 16. The only possible contribution of Mnnm to -/Vfm is estimated to be 
about mc/2M of the A mechanism contributions to Mx and iVf^,. Therefore 
we will neglect the A mechanism for the calculation of in the present 
work. 

We calculate the nuclear matrix elements for Of 2i Ou^f] decay of ''^Ge and 
^°°Mo using the method given in [11]. We describe the initial Oi' and final 2i 
nuclear states in terms of the Hartree-Fock-Bogoliubov type wave functions 
which are obtained by variation after particle-number and angular-momentum 
projection [11-13]. For the case of ^^Ge decay, the calculation of the matrix 
elements Mj with 1 < i < 13 has been performed in [11]. In the present work 
we calculate only the new ones with 14 < i < 25 using the nuclear wave 
functions obtained in [11]. In order to calculate all Mj with 1 < i < 25 for 
the -'^'^°Mo decay, the nuclear wave functions of ^^^Mo{Of) and ^^^Ru{2f) are 
constructed in the same manner as in the case of the ^^Ge decay. Table 2 shows 
the calculated matrix elements for the ^^Ge and ^°°Mo decays as a sum 
of the products CmjMj. It should be noted that the matrix elements of the 
operators with rank spin part, i.e. Mi, M4, M14 and M15 have the dominant 
contributions to Mm. Table 3 summarizes the calculated matrix elements Ma, 
M^, and M^ for the ^^Ge and ^™Mo decays. 

The differential rate for O"*" — > 2"*" Ot//5/3 decay with the energy of one of the 
emitted electrons ei and the angle between the two electrons $12 can be written 
as 

°" = a(°)(ei) + a(i)(ei)Pi(cos^i2) + a(^)(ei)P2(cos ^12) . (9) 



deid cos 6*12 



Each of the angular correlation coefficients a^''\ei) {k = 0,1,2) can be ex- 
pressed as a sum of the products of an electron phase space factor and a 
second order monomial of Zj± defined in Eq. (2). The exphcit form of a^^\ei), 
which yields (ln2)/2 times the right hand side of Eq. (1) upon integration 
over €1, can be readily obtained by combining the relevant equations in [3]. 
Since the expressions for a^^\ei) and a^'^\ei) are rather complicated, they will 
be given elsewhere. Numerical calculations show that a^^\ei) is dominated 
by a term with the factor — (Zi+)^ -|- Z2+Z2- times a positive function of ei, 
whereas a^^\ei) by a term with the factor 2Zi+Zi_ — {Z2+Y — (^2-)^- For 
later reference we denote these two factors as z^^^ and 2;^^) , respectively. 

Figure 1 shows the single electron spectra dWoi^/d^i = 2a*^°) and the ratios of 
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the angular correlation coefficients a'^^^a'^'^^ and a^"^^ ja^^^ for the three limiting 
cases, (a) (A) 7^ 0, (b) {rf) 7^ and (c) {rriu) 7^ 0. Since the coefficients a'^^^ei) 
depend on the parameters (A), {rf) and {rriy) through the results shown in 
Fig. 1 are independent of nuclear models for the cases (a) and (c). We can also 
easily understand the signs of a^^^ and a^^^ from the relations z^^^ = —{Mx{\)Y 
and ^(2) = 2(AfA(A))^ for the case (a), and z^^^ = -2{M„,{m^) /m^f and 
z^'^^ = —A{Mai{m^)/me)'^ for the case (c). On the other hand for the case (b), 
we obtain = -{M^{'q)f + {M^{'n)y and z^^^ = 2{M^{r])f -2{M'^{r])f , imd 
consequently a cancellation between the contributions of and M'^ occurs 
when these are of comparable magnitudes. This is just the case for the ^^'^Mo 
decay, but not for the '^^Ge decay where is much smaller than so that 
there is no significant difference between the cases (a) and (b) in the angular 
correlation. It should also be noted in Fig. 1 that the single electron spectra 
for all the three cases (a), (b) and (c) have approximately the same shape. 
This is in contrast with the 0+ — > 0+ decays where the spectrum for (A) 7^ 
is very different from those for {rrii,) 7^ or (77) 7^ [2,3]. 

Using the matrix elements in Table 3 and the phase space integrals Fj± cal- 
culated in [3], we can deduce from the experimental data r^^i^^ ~^ 2^) 

> 8.2 X 10^3 yr (90% C.L.) [14] for the ^^Ge decay the constraints on the right- 
handed current couplings and the effective neutrino mass listed in Table 4. As 
for the ^°°Mo decay, the Osaka group has obtained the limit Ty2i0~^ ^ 2^) 

> 1.4 X 10^^ yr (68% C.L.) [15] assuming (A) 7^ 0. Because of the differences 
in the angular correlation as we see from Fig. 1, an analysis of the same raw 
experimental data might yield a half-life limit significantly different from the 
above value especially for the case (77) 7^ 0. However we assume here just the 
same half-life limit also for the cases (rj) 7^ and (m^) 7^ in order to compare 
the resulting constraints with those from the ^^Ge data. 

The limits which can be deduced from the experimental bound T[y2(0'^ — > 
0+) > 5.7 X 10^5 yr (90% C.L.) [16] on the 0+ ^ 0+ decay of ^^Ge using the 
nuclear matrix elements of [17] are |(A)| < 3.8 x 10"''', |(?7)| < 2.2 x 10"^ and 
|(mj,)| < 0.19 eV. Comparing these limits with those of Table 4, we notice 
the considerable difference in the absolute sensitivities between the 0+ — > 0"*" 
and 0+ — > 2+ decays, which refiects the smaller Q-value as well as the higher 
electron partial waves associated with the latter. However, it should be stressed 
here that the relative sensitivities to (m,^) and (r)) are comparable in both 
cases. In other words, (rriiy) = 1 eV would give roughly the same decay rate as 
{rj) = 10~^ in the 0+ — > 2+ as well as in the 0+ — > 0"*" decays. At the same time 
it should also be noted that the 0+ — > 2+ decay is relatively more sensitive to 
(A). 

In summary, we have calculated O"*" 01/(3(3 decay rates taking into ac- 

count the recoil corrections to the nuclear currents. As a result, the expression 
for the decay probability becomes a quadratic form of not only the effective 
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coupling constants (A) and (rj) of the right-handed leptonic currents but also 
the effective neutrino mass {rriu) which would be totally absent without the 
inclusion of the recoil corrections. In other words, the recoil corrections give 
the lowest order contribution to the 0+ — > 2+ Oi^PP decay for the case where 
(A) = (77) = and (mj,) 7^ 0. Furthermore, by the numerical calculation of 
the relevant nuclear matrix elements, we have demonstrated that the relative 
sensitivities of 0"*" 2+ decays to (rriu) and {rj) are comparable to those of 
0"*" —>■ 0"*" decays. 

The author thanks S. Yamaji at RIKEN for his warm hospitahty. 
References 

[I] W.C. Haxton and G.J. Stephenson, Jr., Prog. Part. Nucl. Phys. 12 (1984) 409. 
[2] M. Doi, T. Kotani and E. Takasugi, Prog. Theor. Phys. Suppl. 83 (1985) 1. 

[3] T. Tomoda, Rep. Prog. Phys. 54 (1991) 53. 
[4] R. Mohapatra, Phys. Rev. D34 (1986) 3457. 
[5] J.D. Vergados, Phys. Lett. B184 (1987) 55. 

[6] M. Hirsch, H.V. Klapdor-Kleingrothaus and S.G. Kovalenko, Phys. Rev. D53 
(1996) 1329. 

[7] A. Morales, Nucl. Phys. B (Proc. Suppl.) 77 (1999) 335. 
[8] J. Schechter and J.W.F. Valle, Phys. Rev. D25 (1982) 2951. 
[9] J.F. Nicves, Phys. Lett. 147B (1984) 375. 
[10] E. Takasugi, Phys. Lett. 149B (1984) 372. 

[II] T. Tomoda, Nucl. Phys. A484 (1988) 635. 

[12] T. Tomoda, A. Faessler, K.W. Schmid and F. Griimmer, Nucl. Phys. A452 
(1986) 591. 

[13] K.W. Schmid, F. Griimmer and A. Faessler, Nucl. Phys. A431 (1984) 205. 

[14] B. Maier, Nucl. Phys. B (Proe. Suppl.) 35 (1994) 358. 

[15] N. Kudomi et al, Nucl. Phys. A629 (1998) 527c. 

[16] L. Baudis et al, Phys. Rev. Lett. 83 (1999) 41. 

[17] T. Tomoda and A. Faessler, Phys. Lett. B199 (1987) 475. 



6 



(a) (A) 7^0 (b) (77)7^0 (c) (m,)7^0 
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(e -m c2)/Q_ 
^1 e /S/s 

Fig. 1. Single electron spectrum dl^oi//dei in arbitrary units and the ratios of the 
angular correlation coefficients a^^^ /a^'^^ and a^^^ / a^^^ for the 0^ 2f Ou^P decay 
of ^'"^Mo. They are all plotted against the kinetic energy fraction of one of the two 
emitted electrons, where (5/3/3(0+ — > 2^) = 2.494 MeV. Only one of the three lepton 
number violating parameters is assumed to be nonvanishing for each of the three 
cases: (a) (A) / 0, (b) {rj) ^ and (c) {m„) / 0. 
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Table 1 

The operators Minm and the coefficients Cmi , the latter in units of the electron- 

micloon mass ratio Wf/AI. 



i 






1 


-VSh'Soy2 


-i[Ai/3(5'v/5A) + |] 


2 


h'S2 


g[/^/3(5'v/5'A) - 1] 


3 






4 


h'y2 




5 






14 




1 

3 


15 


hy'2 


-{9w/9Af 


16 


HS2 


-kM9v/9A) - 1] 


17 


hS2y'o 


1 


18 


h[S2 0y[]^^^ 


2 


19 


h[S2'S^y'2f^ 




20 


h[S+^y[](^) 


^(ffv/5A) 


21 




^(w/5a) 


22 




1 
4 


23 


/i[Sl®Y'2](2) 


V3 
4 


24 


/i[S'_(8)r'i](2) 


-jTjC^v/s'a) 


25 


/l[5_®r'2](2) 


-i^(5v/5A) 
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Table 2 

Calculated matrix elements Mm for the ^^Ge and ^°^Mo decays. The entries are the 
values of the products CmiMj and their sum Mm in units of lO^^fm"^. 



i 


76Ge 


looMo 


1 


-0.0229 


-0.0077 


2 


0.0013 


0.0003 


3 


0.0002 


0.0008 


4 


-0.0017 


-0.0011 


5 


-0.0003 


0.0007 


14 


-0.0191 


-0.0227 


15 


-0.0128 


-0.0112 


16 


-0.0033 


-0.0006 


17 


-0.0017 


0.0000 


18 


-0.0028 


0.0019 


19 


-0.0005 


0.0001 


20 


0.0006 


0.0005 


21 


0.0002 


-0.0001 


22 


0.0020 


-0.0000 


23 


0.0020 


-0.0026 


24 


-0.0047 


-0.0001 


25 


0.0012 


-0.0010 


sum 


-0.0624 


-0.0427 



Table 3 

Calculated matrix elements for the 0+ — 2'^ Ov^P decays of ^^Ge and -"^^^Mo in 
units of 10~^fm~^. 





Mr, 


M' 


Mm 


^^Ge 1.81 ^ 


13.37 ^ 


0.18 ^ 


-0.0624 


lo^Mo -6.33 


3.38 


5.17 


-0.0427 



^ Ref. [11]. 
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Table 4 

Constraints on the right-handed current couphngs and the effective neutrino mass. 



76Ge lo^Mo 

I (A) I < 8.9x10-'' < 3.9x10-^ 

\{7])\ < 1.2x10-^ < 4.3x10-^^ 

|(to^)| [eV] < 1.0 X 10"^ < 2.2x10^^ 

^ Assuming the same hmit on ^ the (A) mode. 
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